arXiv:1507.05476vl [nlin.SI] 20Jul2015 


The Sylvester equation and the elliptic Korteweg-de Vries system 


Ying-ying Sun 1 , Da-jun Zhang 1 * Frank W. Nijhoff 2 
1 Department of Mathematics, Shanghai University, Shanghai 200444, P R- China 
2 Department of Applied Mathematics, School of Mathematics, University of Leeds, LS2 9JT, UK 


July 21, 2015 


Abstract 


The elliptic Korteweg-de Vries (KdV) system is a multi-component generalization of the 
lattice potential KdV equation, whose soliton solutions are associated with an elliptic Cauchy 
kernel (i.e., a Cauchy kernel on the torus). In this paper we generalize the class of solutions 
by using a Sylvester type matrix equation and rederiving the system from the associated 
Cauchy matrix. Our starting point is the Sylvester equation in the form of kM + Mk = 
rc T — gK~ 1 rc T K^ 1 where k and K are commutative matrices and obey the matrix relation 
k 2 = K + 3ei I + gK _1 . The obtained elliptic equations, both discrete and continuous, are 
formulated by the scalar function SW which is defined using (fc, K, M, r , c) and constitute 
an infinite size symmetric matrix. Lax pairs for both the discrete and continuous system are 
derived. The explicit solution M of the Sylvester equation and generalized solutions of the 
obtained elliptic equations are presented according to the canonical forms of matrix k. 

Keywords: The Sylvester equation, elliptic KdV systems, Cauchy matrix approach, solu¬ 
tions 


PACS: 02.30. Ik, 02.30.Ks, 05.45.Yv 


1 Introduction 


The following elliptic lattice potential Korteweg-de Vries KdV (elpKdV) system 



(1.1a) 

(1.1b) 

(1.1c) 

(Fid) 


(s — s)(w — w) = [(a + u)s — (b + u)s]s — [(a — «)s — (b — u)s\s , 
(s — s)(tc — w) = [(a — u)s + (b + u)s]s — [(a + u)s + (b — u)s]s, 

/ ItV / ^ VO v o 7~. / 

(a + u ——){a — u - ) = a — R(ss), 


s 


s 



(l.le) 


s 


s 


and the continuous elliptic potential Korteweg-de Vries KdV (epKdV) system 



(1.2a) 


s 


s 


A t =4 A 


-XXX 


6 A~A X + 6A x R(s 2 ) - —(R(s 2 )) x , 


(1.2b) 
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with A = —u + y were first derived in [1] through the direct linearisation approach. Here R(x ) 
is associated with the elliptic curve 

y 2 = R(x) = — + 3ei + gx, (1.3) 

x 

where e\,g £ C are moduli of the elliptic curve. In equation (1.1) we use the conventional 
tilde-hat notations to express shifts w.r.t. discrete variables, e.g., 

U — u(n,m) — Un,mi U — %+l,mj U — f-1 5 ^ — ^n+l,m+l• (1-4) 

The direct linearisation approach starts from an integral over certain set of contours, from 
which an infinite order matrix U is introduced. Then the closed forms of the elements of U 
yield nonlinear lattice equations (as examples, see [1-4]). 

In this paper, we will rederive the above two elliptic systems using the Cauchy matrix 
approach, which was successfully applied in [5, 6] to derive integrable lattice equations and to 
analyse their underlying structures. The direct linearisation approach and the Cauchy matrix 
approach look different w.r.t. their individual procedures, but in fact they are deeply related 
together. The former needs a Cauchy kernel as a key auxiliary role and both approaches make 
use of a same infinite order matrix U of which elements generate scaler nonlinear equations. We 
use the notation S for this infinite order matrix in our paper. 

The Cauchy matrix approach is a pure algebraic procedure and it enables us to obtain 
equations, their explicit soliton solutions and Lax pairs. In the Cauchy matrix approach, the 
Sylvester equation 

AM - MB = C (1.5) 

can be viewed as a starting point [7,8]. The matrix M is a dressed Cauchy matrix (see the 
factorization (3.17)) and is used to introduce r-function [5-7]. 

In the present paper we start from the following Sylvester equation: 

kM + Mk = rc T — gK~ 1 rc T K~ 1 , (1.6) 

where r = {r\,r 2 , ■ • • , r/v) r , c = (ci, C 2 , • ■ ■ , cn) t and k, K £ Cjvxiv obey the matrix relation 

k 2 = K + 3dl + gK^ 1 , kK = Kk, (1.7) 

in which I is the IVth-order unit matrix. Based on the above Sylvester equation, the dispersion 
relations for the elpKdV system is defined by 

(al — k)r = (al + k)r , (bl — k)r = (bl + k)r , (1.8) 

and for the epKdV system by 

r x = kr, r t = 4 fc'V, = c T k, = 4 c T k 3 . (1.9) 

In next section, we will focus on the matrix system (1.7). We will see that the matrix 
system (1.7) is actually governed by the points on the elliptic curve (1.3). Sec.3 will focus on 
the Sylvester equation (1.6) and the scalar function S™' defined in (3.28). Explicit solution 
M of (1.6) will be given for the cases of k being diagonal and Jordan block form and their 
combination. Solutions of the elliptic elpKdV and epKdV equations are consequently obtained. 
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As a generic element, composes an infinite order matrix S (same as the matrix U in the 

direct linearisation approach in [ 1 ]). satisfy some recurrence relations which can be 

viewed as discrete equations of Shd) defined on Z x Z and will play crucial roles in deriving the 
continuous epKdV system. 

Apart from the Sec .2 and 3 above mentioned, in Sec.4 and Sec.5 we respectively derive the 
elpKdV system and epKdV system together with their Lax pairs. In Sec .6 we discuss continuum 
limits of the elpKdV system. Sec. 7 is for conclusions. Besides, in Appendix we list properties 
of lower triangular Toeplitz matrices which play important roles in our paper. 


2 Points on the elliptic curve: parametrization and selection 

2.1 Scalar case 


Consider the elliptic curve 

fe2 = jR (^) =K + 3ei + f- 

The discrete plane wave factor is defined as 


Pi = 


f a + kj \ n 
\a-kij 


b + h \ m o 

b^kj 


( 2 . 1 ) 


( 2 . 2 ) 


with a phase factor p®, where we require that ki together with Ki obeys the elliptic curve ( 2 . 1 ), 

i.e. 

k$ = Ki + 3ei + J-. (2.3) 

■K-i 

In the case of classical soliton solutions (see [6,9]), {Lj} play the role of wave numbers, which 
should be distinct so that they can represent different solitons. Since ki and Ki are coupled 
through the elliptic curve (2.3), (say, both ki and —ki correspond to a same Ki ,) we consequently 
require that they can identify each other, i.e. 


h / kj Ki / Kj. (2.4) 

Note that for the arbitrary two points ( ki,Ki ) and ( kj,Kj ) on the elliptic curve (2.1) we always 
have the relation 

(ki + kj)(ki - kj) = (Ki - Kj) Ki *y~ 9 . (2.5) 

This means that if we take 

(ki + kj)(KiKj - g) / 0, (2.6) 

then (2.4) is guaranteed. Equation (2.6) is the criteria that we select points from the elliptic 
curve ( 2 . 1 ). As a consequence of ( 2 . 6 ), ki / 0 . 

The elliptic curve (2.1) can be parameterized using Weierstrass’s elliptic function p(n) as 
the following (cf. [ 1 ]): 

K = p(n)-e r, k= f[ K) (2.7a) 

2 (p (k) - ei) 

ei = p(w), 5 = (ei-e 2 )(ei - e 3 ), (2.7b) 
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where ei = p(w + u/), e 3 = p(u/) and w and u/ are respectively the half periods of p(k). 

An elliptic function is a meromorphic function with double (complex) periods. Suppose that 
f(n) is an elliptic function and its double periods are 2u and 2u/, respectively. Let B denote a 
fundamental period parallelogram ABCD as described in Fig.l. The border segments AB and 
BC together with points {A, B,C} are not included in B due to the double periodicity. The 

—id id id 1 id -f- id 



Figure 1: Fundamental period parallelogram B 

elliptic function /(re) obeys the Liouville’s theorems (cf. [10]) that is, for the elliptic function 
/(re), in the period parallelogram B the times that it takes value a equals to the order of /(re). 
Under the parametrization (2.7), for the points on the curve (2.1) we have 

Ki = p(^i) — ei, k i= f /(ACi) (2.8) 

2 (p(«i) - ei) 

Then the criteria (2.6) can alternatively be described through the following requirement on Kf. 



(2.9a) 

(p(«t) - e i)(p(nj) ~ ei) ^ g, 

(2.9b) 

“h fcj 0 ? 

(2.9c) 


for i,j = 1,2, • • • ,N. We note that here and after when we talk about Ki + Kj we always mean 
Ki + Kj modulo the periodic lattice, i.e. Ki + Kj{ mod(2w, 2u/)). In fact, 0 is the pole of p(re) and 
thus we require Ki 7 ^ 0 to avoid singularities. Besides, ki can not be zero as a consequence of 
ki + kj 7 ^ 0. In the light of the Liouville’s theorems (cf. [10]) since p'(k) is a third-order elliptic 
function in the period parallelogram B it only has 3 zeros which are and 00 + u/. To avoid 
breaking the one-to-one correspondence of ki and Ki, we require Ki ^ + u/}. Thus we 

have (2.9a). (2.9b) is from the requirement KjKj 7 ^ g. For the third one in (2.9), we can prove 
that under (2.9a) and (2.9b) we have 

Ki + Kj = 0 ki + kj = 0 . ( 2 - 10 ) 

In fact, since p(re) is even and p'(k) is odd, from the parametrization (2.8) we immediately find 
that if Ki + Kj = 0 then ki + kj = 0. On the other hand, under (2.9b), from the factorization 
(2.5), if ki + kj = 0 and (2.9b) holds, there must be Kj = Kj, which means Ki = ±Kj in B in the 
light of Liouville’s theorems (cf. [10]). The case Ki = Kj is impossible because this case yields 
ki = kj = 0 due to ki + kj = 0. ki = 0 requires p'(Ki) = 0 which is impossible in B'. Thus, 
Ki = — Kj is the only choice, i.e. Ki + Kj = 0. 


4 






2.2 Matrix system (1.7) 

Let us come to the matrix relation (1.7). Although at first glance (1.7) suggests an interpre¬ 
tation of this matrix relation as a matrix version of an elliptic curve, it just represents the 
coordinatization of a collection of points on the given elliptic curve ( 2 . 1 ). 

To understand this, let us consider a similarity transformation 

ki = TkT 1 , Ki = TKT 1 , (2.11) 


where T serves as the transform matrix. Obviously, under the above similarity transformation, 
(1.7) is formally invariant: 

fcy = K\ + 3eiJ + gKi \ k\K\ = K\k\. (2.12) 

Thus, in the following we only need to consider the relation 

T 2 = K + 3e 1 I + gR- 1 , TK = KT (2.13) 

where T is the canonical form of k. 

When 

T = Diag(fci,fc 2 ,--- ,k N ), (2.14a) 

K is taken as 

K = Diag(A 1 ,i7 2 , • • • , K n ), (2.14b) 

where {ki,K{) are the points on (2.1), i.e. satisfying (2.3). {ki} is the eigenvalue set of k. Here 
we request that each fc, / 0 and kf ^ k'j for % ^ j. Under such a requirement one can see that 
the criteria (2.6) is satisfied in the light of the factorization (2.5). 

It is interesting to consider the case that T is a IV-th order Jordan block 


r = 


/ h 


l 

o 


0 0 
k\ 0 
1 k\ 


0 0 \ 
0 0 
0 0 


k\ 7 ^ 0 . 


\ 0 0 0 ••• 1 ki / 


(2.15) 


In this case, k\ is the only eigenvalue of k with algebraic multiplicity N (and geometric multi¬ 
plicity 1). In terms of the parametrization (2.7) we need k\ satisfies 


ki <E B\ {0,w,w',a; + w'}i p(«i) / V5 + ei- 


(2.16) 


To find K that corresponds to the Jordan block (2.15), we will make use of properties of Toeplitz 
matrices. Let us shortly introduce such matrices. For more details please see Appendix A. A 
matrix of the following form 




a\ 

0 

0 

• 0 

0 

\ 


C12 

o\ 

0 

• 0 

0 



0,3 

02 

o\ 

• 0 

0 


(2.17) 

ON 

ON-1 

ON-2 • 

• a 2 

a\ 

/ 
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is called a lower triangular Toeplitz (LTT) matrix. All the IV-th order LTT matrices compose a 
commutative set T [Jv h The complex value matrix (2.17) can be generated by an analytic function 
f(k) at certain point k = k 0 through 

q3~ 1 

a j = jj^Yy f( k )\k=k 0 , (i = 1,2,--- ,AT). (2.18) 

With this correspondence we denote (2.17) by T [N] [f(ko)], and we have properties (See proposi¬ 
tion A.2) 

(TW[f(k 0 )]) 2 = TW[f(k 0 )}, = T [N] [l/f (ko)\. 

We also note that the Jordan block (2.15) is a LLT matrix generated by f(k) = k at k = k\. 

Now we can look for a LTT matrix K which satisfies (2.13) with T (2.15). In fact, according 
to proposition A.l, for the Jordan block (2.15), when TK = KT there must be K € T [A L 
By taking derivatives with respect to k of the elliptic curve (2.1) at the point (k±,Ki) where 
K\ = K(k\ ) and K(k ) is viewed as an implicit function of k determined by the curve (2.1), we 
find (for i > j) 

= + + i^k = k 1 ). (2.19) 

The l.h.s. and the first term and third term on the r.h.s. respectively correspond to the elements 
of the LTT matrices generated by f(k) = k 2 , K(k) and 1/K(k) at k = k\. This means for the 
Jordan block case (2.15) if we take K = T [Ni [K(ki)\ then the relation (2.13) holds. 

In Sec.3.2.3 we will list solutions of (2.13) using unified notations. 


3 The Sylvester equation and infinite matrix structure 

In this section we will first investigate solutions of the Sylvester equation (1.6) and derive explicit 
expression of M. Then, with the help of some special matrices we will investigate recurrence 
relations of scalar function (defined in (3.28)) and properties of the infinite matrix S 

composed by S^\ 

3.1 Solvability of (1.6) 

For the solution of the Sylvester equation (1.5), there is the following well known result [11]. 

Proposition 3.1. Denote the eigenvalue sets of A and B by 6(A) and 6(5), respectively. For 
the known A,B and C, the Sylvester equation (1.5) has a unique solution M if and only if 
6 (A) n 6(B) = 0 . 

Based on this proposition, we find the following. 

Proposition 3.2. Consider the Sylvester equation (1.6), i.e. 

kM + Mk = rc T -gK~ 1 rc T K- 1 , (3.1) 

where the matrices k and K satisfy 

6(fc)n6(-fe) = 0, 

Z(gK~ l )cZ{K) = 0 , 
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(3.2a) 

(3.2b) 




and the matrix relation (1.7), i.e. 


k 2 = K + 3eiJ + gK- 1 , kK = Kk. 
Then, the ‘dual’ matrix equation 


KM - MK = k rc T — rc T k 


(3.3) 


(3.4) 


holds. 

Proof. Here we note that the condition (3.2a) is necessary to guarantee the solvability of equation 
(3.1) in light of proposition 3.1. Then, left multiplying k to (3.1) yields 

k 2 M + kMk = k(rc T — gK~ 1 rc T K~ 1 ). 

Next, using the Sylvester equation (3.1) and replacing the term kM with —Mk + rc T — 


gK- 1 rc T K~ _1 , one has 

k 2 M - Mk 2 = —rc T k + krc T + gK^r^K^k - gkK^r^R- 1 . (3.5) 

Making use of the relation (3.3) to replace k 2 we find 

gK- 1 (KM - MK - krc T + rc T k)K~ l = KM - MK - krc T + rc T k. (3.6) 

This can be rewritten as a Sylvester equation 

gK 1 W-WK = 0, W = KM - MK - krc T + rc T k. (3.7) 

Based on proposition 3.1 and noting that £(gK~ 1 ) n £(K) = 0 in (3.2b), the equation (3.7) 
has a unique solution W = 0, which means (3.4) holds. □ 


Here we note that the condition (3.2) is natural because it is actually the criteria (2.6) for 
selecting points from the elliptic curve (2.1). We also note that we cannot derive (3.1) from 

(3.4) . In fact, we start from (3.4), replace K using (3.3) and we get 

k 2 M - Mk 2 - g(K~ l M - MK~ 2 ) = krc T - rc T k. 

At meantime, from (3.4) we also have 

K~ 2 M - MK 1 = —K~ 1 krc T K~ 1 + K~ l rc T kK~ l . 

Combining them together and using the commutative relation kK = Kk we reach 
kY-Yk = 0, Y = kM + Mk-rc T + gK 1 rc T K~ l . 

Obviously, Y = 0 is a solution to the above equation but it is not unique. This means (3.1) and 

(3.4) are not equivalent. (3.1) is more general and (3.4) is a by-product of the former. In the 
following discussion it is sufficient that we only consider (3.1). 
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3.2 Solution to the Sylvester equation (3.1) 

3.2.1 Canonical form of (3.1) 

Using the similarity transformation (2.11) and denoting 

Mi=TMT ~\ n=Tr, cl = c T T~ 1 , 

it follows from (3.1) that 

M\k\ + k\M\ = r\ c[ — gK^ ] r\c(K^ 1 , 

which is the same form as (3.1). This means when we solve the Sylvester equation (3.1) 
need to consider the following canonical form 

MY + YM = rc T - gK^r^K -1 , 

together with (2.13), i.e. 

Y 2 = K+ 3 ei I + gR- 1 , YK = KY, 

where 

r = (n,r 2 , ■ ■ ■ ,r N ) T , c = (ci,c 2 , • • • ,c N ) T 
and we suppose Y is the canonical form of k. 


3.2.2 List of notations 

Based on Sec.3.2.1 we only need to consider three basic cases of T: 


Cd ({fcjli ) = Diag(A: 1 , fc 2 , ■ 
/ k\ 0 O' 
1 ki 0 ■ 

0 1 h ■ 


Y^ih) = 


,k N ), ( k 2 ^ k 2 , ki / 0), 

0 0 \ 

0 0 

' =T M N, 


0 0 


0 0 0 ••• Ik! ) 

r™ = Diag^^d^'lf 1 ), r^^fcjVi+i), rf 3] (A: J v 1 + 2 ), ■ ■ ■ , r^ ] (^ 1+ (-i)))’ 


(3.8) 

we only 
(3.9a) 

(3.9b) 

(3.9c) 

(3.10a) 

(3.10b) 

(3.10c) 


where Ylj=i^j = N. The subscripts d, j and q correspond to the cases of Y being diagonal, 
being of Jordan block and generic canonical form, respectively. For convenience of the later 




discussions, let us collect some notations below. 


JV-th order vector : e [N] = (1,1,1, • • • , 1) T , 
JV-th order vector : e [ ,' vl = (1,0,0, • • • , 0) T , 


IV-th order vector : g^(a) = 


-1 1 

9~ ) q" 5 * * 

a z a^ 

(-1) 

’ a 7 

N x N matrix : 


= ( Gi,j)NxN, 

Gi,j = 




f Cl ■■■ 

Ci 

1 

to 




c 2 • • • 

Cat-1 

N x N matrix : 


= 

c 3 ■ ■ ■ 

CN 




\ CN ■■■ 

0 


1 - g/jKjKj) 
k^ “I - kj 


CN 0 

0 0 

0 0 / 


,}j v ) is called a skew LTT, which is also introduced in Appendix A. 


(3.11a) 

(3.11b) 

(3.11c) 

(3.lid) 

(3.lie) 


3.2.3 Solutions to the matrix system (2.13) (i.e. (3.9b)) 

We have discussed solutions to the matrix system (2.13) (i.e. 3.9b) in Sec.2.2. Here we list them 
using the notations given in Sec.3.2.2. 

Proposition 3.3. The matrix system (2.13) admits the following three cases of solutions: 


(1) . Diagonal case: 

r = r^({%}f), k = ( 3 . 12 ) 

where 

k] = Kj + 3ei + gK~\ j = 1,2, • • • , N. (3.13) 

(2) . Jordan block case: 

r = rf 7 (fci), K = T lN, [K{k l )\. (3.14) 

(3) . Generic case: 

r = r 7 G 7 , (3.15a) 

K = Diag(r^ 7 ({^}f 1 ),T^ 7 [^(fc JVl+ i)],--- ,T^[K(k Nl+{a _ 1) ])). (3.15b) 


3.2.4 Solutions to (3.9a) 

Now let us come to the solutions to the Sylvester equation (3.9a). 

Case i. r = r^ ] ({%} 7V ). 

Solution to (3.9a) is given by 

M - FG^lk^H - 

where 

F = Diag(ri,r 2 , • • • ,r N ), H = Diag(ci,c 2 , • • • ,c N ) 


(3.16a) 

(3.16b) 
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Case 2. T = T l J N] (k 1 ). 

This is also referred to as the Jordan block case. In this case, T and K take (3.14). We 
suppose 

M = FGH. (3.17) 

where 

F = T^({r,-}f), H = HW({ Cj }?), (3.18) 

and G is a N x N unknown matrix, r and c can be expressed through F and H as 


r = F e i, c = H e\, 

where ei = is defined in (3.11b). Then, the equation (3.9a) is written as 
FGHT + TFGH = Fei el H - gK^ 1 F e x ef HK 1 . 
Further, thanks to proposition A.4, one has 


FGY T H + FTGH = F ei e[ H - gFK 1 ei e[ K~ lT H, 


and then 

To solve it, we set 


GT T +TG = e l el-gK~ l e 1 elIC 


-1 T 


G = (Gi,G 2 ,--. ,G n ) 

with column vectors {Gj}. (3.22) is then written into the following equation set, 

(fc 1 / + r)Gi=e 1 --|-A 1 , 

(hi + T)Gj + \ + Gj = -|(< ^)Ai, (j = 1,2, ■ ■ ■ , N - 1), 


where K\ = I\(h) and 

Ai = (ai, a 2 , • • • , oat) t , aj = — 
The above is solved by 


1 - 7 -l 1 


:di~ 


J 0-1)! AV 


(j = 1,2,- -- ,N ). 


Gj = 


di l g [N] (a) |a=2fci 


l~ l + g E(yr^(C‘jO r "( 2t '> ' Au 0 = 1 - 2 '"' < N >- (3.25) 


(3.19) 

(3.20) 

(3.21) 

(3.22) 

(3.23) 

(3.24a) 

(3.24b) 

(3.24c) 


Cases. r = r^. 

In this case, we still suppose the factorization (3.17), where 
F = Diag(r ^ l] ({kj ), T [iV2] ({rj ), • • • ,T™({r J -}^+^+;;^_ 1+1 )), (3.26a) 

fr = Diag(ff™({^}fO,«f al ({c J -}Kf 2 ),--- ,^f s] ({c,'}SSt:S_ 1+1 ))i (3.26b) 

G is a symmetric matrix with block structure 

G = G T = (Gijsxs (3.26c) 
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and each Gjj is a JVj x Nj matrix. It is not difficult to find G t j(= Gjf) are given by 


= }f), 

Gjj = (Gn,Gi 2 , ■ ■ ■ , GiNj), (1 < i < j < s), 


(3.27a) 

(3.27b) 


with 




(- 1 ) 


l—i 


^ (l-i)\\ dk K(k) 


l—i 




k=k 


JVi+l 


[«i] 


Gil = 


( 1 - 1 ) 

l 


(- 1 ) 


m , 1 

i l—m 1 


_ 

(Z — m)! V k K(k) 




/c fciv^+(j_l) 


for 1 < i < j < s, where a, = fcj . +fcjv | +0 _ 1) , fiij = k Nl+ ^_ X ) + £^+ 0 - 1 ) and k and K{k) satisfy 
the elliptic curve (2.1), i.e., k 2 = K 2 (k) + 3ei + g/K(k). 


3.3 Infinite matrix S 

3.3.1 Recurrence relation of S^ 1 ’^ 

Let us go back to the Sylvester equation (3.1) and the matrix relation (3.3). Using the elements 
{M,k,K,r,c} in (3.1) and (3.3), we introduce an oo x oo matrix S — (*S^ , ^)ooxoo; LJ £ 1^, 
where the elements are defined as(cf. [1]) 


s ( 2 i, 2 j) = c t K 3 ^j + M)- x K L r , (3.28a) 

s ( 2 i+i, 2 j) = c T K j(i + M)~ l kK l r, (3.28b) 

S ( 2 i, 2 j+ 1 ) = c T K j k (j + M)~ 1 K i r, (3.28c) 

s ( 2 i+i, 2 j+i) = c T K j k (j + M)~ 1 kK i r. (3.28d) 


For these elements we have the following. 

Proposition 3.4. For the saclar functions defined in (3.28) with (M, K , k, r, c) satisfying 
the Sylvester equation (3.1) and the matrix relation (3.3), we have the following relations, 


S— 1 

g{hj+2s) — _ 'y _ g(2s—2l—2j)g(i,2l+l)^ 

1=0 

s 

g(i,j-2s) _ g(i-2s,j) _|_ ^^(g(-2s+2l-l,j)g(i,—21) _ g(2l—2s-2,j)g(i,-2l+l)^ 
1=1 


(3.29a) 

(3.29b) 


where s = 1, 2, • • • . In particular, when s = 1, one has 

s (i,j+2) = s {i+2,j) _ s (i,j)g(i,o) + s (o,j) s (i,i) , (3.30a) 

5 (*,i-2) = s (i-2,j) + s (-ij) s (i,-2) _ s (-2,j) s (i-1 )_ (3.30b) 
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Proof. First, from the Sylvester equation (3.1) we have the following relation 

S — 1 

k s M - (-1 ) s Mk s = J2(- 1 ) j k s ~ j ~ 1 (rc T - gK- 1 rc T K~ 1 )kP , (s = 1,2, •• • )• (3.31) 

3=0 

In fact, the Sylvester equation (3.1) itself is the case of s = 1 of (3.31), while (3.5) is the case 
of s = 2. Making use of mathematical inductive approach we can reach (3.31). Similarly, from 
(3.4) one has a parallel result 

S — 1 

K S M - MK S = Y K s ~ j - l (krc T - rc T k)K j , 

3 =0 
or 

s 

MK S - K S M = Y K~( s ~ j+1 \krc T - rc T k)K~ j , (a = 1,2, ■ ■ ■ )■ 

3 = 1 

Now let us prove the relation (3.29a). We introduce the auxiliary vectors 

u ( 2< > = (J + M)~ l K l r, u (2i+1) = (I + M)“ 1 fcKV, i € Z. 

From this we immediately have 

K s u {2i) + K s Mu {2i) = K s+i r , 

A' s u (2i+1 ) + rMn (2i+1) = kK s+i r. 

Replacing K S M using the relation (3.32a), one hnds 

S— 1 

(I + M)K s u^ = K s+i r - Y K s ~ l ~ l {krc T - rc T k)K l u^ 2i \ 

1=0 

s —1 

(J + M)K s M (2i+1) = fcK s+ V - ^ - rc T fc)K 2 u (2i+1) . (3.35b) 

1=0 

These relations, left-multiplied by c T K^{I + AT) -1 , yield 


(3.32a) 

(3.32b) 

(3.33) 

(3.34a) 

(3.34b) 

(3.35a) 


S — 1 


g(2i,2j-\-2s) ^(22+25,2^') ^ ^ ^^(2s—21—l,2j) £j(2i,2l) ^»(2s—2Z—2,2j) ^(2z,2/H-l) ^ 

Z=0 

s—1 

g(2z+l,2j+2s) _ g(2z+2s+l,2j) _ (2s-2Z-l,2j)g(2t+l,2Q _ g(2s-2Z-2,2j)g(2z+l,2Z+l) ^ 

Z=0 

and left-multiplied by c r IT J fc(/ + AT) -1 , yield 


a—1 

g(2i,2j+2a+l) _ g(2i+2s,2j+l) _ y~^ (g (2a-2i-l,2j+l)g(2i,2i) _ ^(2s-2i-2,2j+l)g(2i,2Z+l) ) 

1=0 

s —1 

j g(2i+l,2j+2s-|-l) _ g(2j+2s+l,2j+l) _ yyg(2s-2l-l,2j+l)g(2i+l,2l) _ g(2s-2l-2,2j+l) g(2i+l,2l+l)^ 

1=0 
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They are merged to the relation (3.29a). 

The relation (3.29b) can be proved in a similar procedure, in which we use the following 
counterpart of (3.35): 

S 

(I + M)K- s u ( - 2l) = K~ s+i r + Y iT~ (s “ /+1) ( krc T - rc T k)K~ l u ^ , (3.36a) 

1=1 

S 

(I + M)K- S u^ 2i+1) = kK s+i r + Y K-( s ~ l+1 \krc T - rc T k)K~ l u^ 2i+1 \ (3.36b) 

i=i 

with s = 1,2,---. We note that (3.30) are corresponding to the algebraic relations (2.13) in 
Ref. [1], □ 

3.3.2 Invariance and symmetry property of S'bd) 

The matrix S (or the element S'*’") keeps invariant under the similarity transformation (2.11). 

We have shown that the matrix relation (1.7) and the Sylvester equation (3.1) keep invariant 
formally in terms of the similarity transformation (2.11) and notations (3.8). Using (2.11) and 
(3.8) one can rewrite (3.28) and find 


S ( 2i ’ 2j) = c[K{(I + Mi)- 1 Kin, 

(3.37a) 

S ( 2 i+I, 2 j) = c[K{(I + M 1 )~ 1 k 1 Kfr 1 , 

(3.37b) 

S ( 2 i, 2 j+ 1 ) = c T K J ki (j + Mi)- 1 Kjn, 

(3.37c) 

s (2i+i,2j+i) = + AI 1 )~ 1 k 1 Kir 1 . 

(3.37d) 


In addition, we have the following symmetry property. 

Proposition 3.5. Suppose that M , K , k,r, s satisfy the Sylvester equation (3.1) together with 
the matrix system (3.3) in which £(fc)n£(— k) = 0 and &(gK~ l ) n £(K) = 0. Then the scalar 
elements defined by (3.28) satisfy the symmetry property 

S (hj) = sM, (3.38) 


i.e. the infinite matrix S is symmetric. 

This can be proved using the property given in proposition A.4 and following the procedure 
described in Appendix C in Ref. [7]. Here we skip the proof. 

Hereafter we always require that {M,K,k,r,s} satisfy the assumption of proposition 3.5, 
under which we proceed further discussions. 

4 The elliptic lattice potential KdV system 

4.1 Discrete dispersion relation and recurrence relations 

Now let us insert discrete dispersion relation on r as follows 

(al — k)r = (al + k)r, (bl — k)r = (bl + k)r, a,6 ^£(±fc), (4.1) 

while we take c to be a constant vector. 
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Employing a similar procedure as done in [8], from the Sylvester equation (3.1), matrix 
system (3.3) and the dispersion relation (4.1), one can derive the following shift relation of M , 


(al — k)M = (aI + k)M , 

(4.2a) 

(bl — k)M = (bI + k)M, 

(4.2b) 

M(al + k)~ (al + k)M = rc T - gK- 1 rc T K~ x , 

(4.3a) 

(al - k)M - M(al - k) = rc T - gK- 1 rc T K~ 1 , 

(4.3b) 

M(bl + k)~ (hi + k)M = rc T - gK^f^R- 1 , 

(4.3c) 

(i bl - k)M - M(bl -k) = rc T - gK- 1 rc T K~ 1 . 

(4.3d) 


These relations lead to the following results. 

Proposition 4 . 1 . Under the assumption of proposition 3.5 and the dispersion relation ( 4 . 1 ), 
the scalar functions defined by ( 3 . 28 ) enjoy the following recurrence relations: 

a g(2i,2j) _ g(2i,2j+l) _ a g(2i,2j) _|_ g(2i+l,2j) _ 5(250) 5(0,25) q. gS( 2 i,-2) g(-2,2j) ^ ( 4 . 4& ) 
a g{2i,2j) _|_ 5(2525+1) _ a g(2i,2j) _ 5(24+1,25) _|_ 5(250) 5(0,25) _ ^g(2i,-2) fj(-2,2j) ^ (^^ 
bS( 2 i ’ 2 j) — _ ^c<( 2 i, 2 j) _|_ 5(24+1,25) _ 5(250)5(0,25) _|_ ^^(2i,-2)^(-2,21) ^ ( 44c ) 

55(2421) _|_ 5(2521+1) _ 55(2521) _ 5(21+1,21) _|_ 5(250)5(0,21) _ ^5(25-2)51-2,21)^ (4 44) 
a g(2i+l,2j) _ 5(21+1,25+1) _ a g( 2 i+l, 2 j) _|_ 5(21+2,21) _ 5(21+1,0)5(0,21) 

+ gS( 2i ~ 2 ’ 2j) + <?S (2i+1 ’-2)51—2,21) + 3ei 5(2521) ^ (4.4 e ) 

a 5(2i+l,21) _|_ 5(21+1,21+1) _ a 5(2i+l,21) _ 5(21+2,21) _|_ 5(21+1,0)5(0,21) 

- gS( 2i ~ 2 ’ 2j) - 5 5(2i+! ~ 2 ) S {~ 2 , 2 j) _ 3^5(2521) ^ ( 4i4 f) 

55(2i+l,21) _ 5(24+1,25+1) _ 55(21+1,21) _|_ 5(24+2,21) _ 5(24+1,0)5(0,21) 

+ gS i2i ~ 2 ’ 2j) + 5 5 ( 2i+! - 2 ) 5(—2,21) + 3ei 5(252!) ? (4.4g) 

55(24+1,21) _|_ 5(24+1,21+1) _ 55(24+1,21) _ 5(24+2,21) _|_ 5(24+1,0)5(0,21) 

- g s( 2i ~ 2 ’ 2j) - gS {2i+1 - 2 ) 5 (- 2 , 21 ) _ 3 ei 5 ( 2421 ) _ ( 4 . 4 h) 

The proof of this proposition is similar to the one for theorem 2 in Ref. [8]. Here we skip 
the details. We also note that these relations are corresponding to the discrete matrix Riccati 
type of relations ( 2 . 12 ) in Ref. [ 1 ]. 

4.2 Elliptic lattice equations 

To obtain elliptic lattice equations, we introduce scalar functions (cf. [ 1 ]) 

u = S<™\ s = S^~ 2,0 \ h = St~ 2 ’- 2 \ v = l iu = l + S(- 2 ’ 1 ). ( 4 . 5 ) 
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It then follows from (4.4) that 

a(u — u) = — (iS^ 0,1 ) + S'*' 0,1 -’) — gss + uu, (4.6a) 

a{h -h) = “ 2) + 5( -1 > _2) ) - ghh + ss, (4.6b) 

a(s — s) = ghs + w — v — us, (4.6c) 

a(s — s) = ghs + w — v — us, (4.6d) 

a(v — v) = S^ 1 ’ 1 ) + 3eis + uv + gvh + gs(S^~ 2 ’~ 1 " 1 + S’*- -2 ’ -1 ^), (4.6e) 

a(y — v) = + 3eis + uv + gvh + g’s(S^~ 2 ~ 1 ' > + S^~ 2, ~^), (4.6f) 

a(w — w) = S^ l) + 3ei s + uw + ghui — s(S^ 0,1 ^ + S^ 0,1 )), (4.6g) 

a(w — w) = + 3eis + uw + ghw — s(5^ 0,1 ^ + 5^ 0,1 ^), (4.6h) 


where we have made use of the symmetric property = S 1 -- 5 ’®) and the relation (3.30b) with 
(i,j) = (0,-2), i.e. 

^(0,-4) = 5 (-2,-2) + 5 (-l,-2) 5 (0,-2) _ £(-2,-2)£(0,-1) = ^ + S|S -(-2-1). 

One more relation coming from (3.30b) is 

sS 1 ! -1 ’ 1 ) = 1 — vw, (4.6i) 

when taking (i,j) = (0,1). The shift relations w.r.t. { b ,T } can be obtained by interchanging 
the relations {a,~} with {5,~ } in (4.6). 

With these relations in hand, we can combine them, eliminate S'! 0 ’ 1 ), 2 ), h and v and 

reach the closed form of elliptic lattice equations of the variables u, s and w, (cf. [1]) 

(a + b + u — u){a — b + u — u) = a 2 — b 2 + g(s — s)(s — s), (4.7a) 

(s — s)(w — w) = [(a + u)s— (b + u)s]s— [(a — u)s — (b — u)s]s, (4.7b) 

(s — 's)(u) — w) = [(a — u)s + (b + u)'s\s— [(a + u)s + (b — u)s]s, (4.7c) 

(a + u — ^)(a — u + —) = a 2 — R(ss), (4.7d) 

s s 

111 111 

(b + u — —){b — u-\ -) = b 2 — R(ss), (4.7e) 

s s 

where 

y 2 = R(x) = — + 3e + gx. (4.8) 

x 

This is an elliptic generalization of the lpKdV equation, namely that if g = 0 the first equation 
(4.7a) decouples and becomes the standard lattice potential KdV equation. 

Furthermore, solutions to the elpKdV system (4.7) can be expressed in the explicit structure: 


5 (°,°) = c T {I + M)- l r , 

(4.9a) 

5 (- 2 ,o) = c t (j + M)~ l K~ 1 r, 

(4.9b) 

g(- 2 - 2 ) = c T K~ x {I + M)~ l K~ l r, 

(4.9c) 

1 _ 5 ( " 1 ’ 0) = 1 - c T {I + M)~ l kK~ l r , 

(4.9d) 

1 + = c T k(I + M)K~ 1 r. 

(4.9e) 
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We only need to solve out r from (4.1) with k taking three cases in (3.10). We will list out 
explicit forms of r in Appendix B. 


4.3 Lax pair for the elpKdV system (4.9) 

Rewriting (3.33) as 

Kr = (I + M)u( 2i \ (4.10a) 

kK l r = (I + M)u (2i+1) , i € Z, (4.10b) 

we can perform a tilde shift to (4.10a) and multiply the result by (al — k) to give: 

K\al + k)r = (J + M)(al - k)u (2i) + (rc T - ff jr 1 rc T JT 1 )u^ ) (4.11) 

where we have also make use of the fact kK = Kk, the dispersion relation (4.1), shift relation 
(4.2a) and the Sylvester equation (3.1). This further brings 

(al - k)u^ 2i) = a« (2i) + tt (2i+1) - S (2i ’ 0) u (0) + gS (2i ’~ 2) u ( ~ 2) . (4.12) 

Taking i = 0 it reads 

(al — k)u^ = (a — u)u ^ + u ^ + gsK~ 1 (wu^ — su ^), (4-13) 

where we have replaced u^ 2 ' 1 by 

u ( - 2) = K~ l (wu (0) - sm (1) ) (4.14) 

coming from the relation (3.35a) with s = 1, i = — 1. 

In a similar way, from (4.10b) we have 


(■ al—k)u ( - 1) = ( K+gwwK 1 )u^+(3ei+gss+a(u—u)—uu)u^+(a+u—gwsK 1 )i^ 1 \ (4.15) 


where 

u® = Ku® - +mx (1) (4.16) 


is used. 

It is easy to check that under the similarity transformation (2.11), equations (4.13) and (4.15) 
are formally invariant if we define u\ ' 1 = Tu^°\ tt^ 1 * = Tu^\ That means we can directly 
consider that k is in its canonical form and K is consequently defined by (3.3). Thus, the first 
row of k will be (k, 0, 0, • • • ,0) and the first rows of K and K 1 have to be (if, 0,0, ■ ■ ■ ,0) and 
(1/AT, 0,0, • • • ,0) respectively, where (k, K) obey the elliptic curve (2.1). 

Let us denote the first element of u^ by (i/°))i and the first element of id 1 ) by (wd)) 1; and 
introduce the vector 




(4.17) 


Then, from (4.13) and (4.15) and their (b,^) version we obtain the following discrete linear 
system: 


(a — k)(j) = L(K)<f>, (4.18a) 

(b — k) 4 > = M(K)(j), (4.18b) 
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where 


(4.19) 


L(K) 


a — u + jfSiv 

K + 3e + gss + a{u — u) — uu + -j^ww 


1 ^ 

a + u — j^ws J 


which is as same as in [1], and M(K) is the (6, counterpart of L(K). The point (k, K ) obeys 
the elliptic curve (2.1) and here they play the roles of spectral parameters. The compatibility 
condition 


LM = ML 


(4.20) 


yields the whole ellpKdV system (4.7) with the exception of eq.(4.7c). Here we nota that there 
can be mis-matched between a discrete system and the system derived from its Lax pair, for 
example, the discrete Boussinesq-type equations [4], 


5 The elliptic potential KdV system 

In this section we try to derive a continuous elliptic potential KdV (eqKdV) system. The 
procedure will be similar to the one for the KdV system in [7], which can be viewed as a 
continuous version of Cauchy matrix approach, where the recurrence relation of ( see (2.3) 
in [7] and (3.30) in this paper) will play key roles. 

5.1 Evolution of M 

We suppose that M, r, c are functions of (x, t) while k is still a non-trivial constant matrix. 
The dispersion relation is now defined through the evolution of r and c as follows, 

r x = kr, c x = k T c, (5.1a) 

r t = 4 k 5 r, c t = 4(fc T ) 3 c. (5.1b) 

Taking the derivative of the Sylvester equation (3.1) w.r.t. x and making using of (5.1a) we 
have 

kM x + M x k = r x c T + rc x — gK~ l r x c T K~ 1 — gK~ 1 rc x K~ 1 

= krc T + rc T k — gK l krc T K 1 — gK~ 1 rc T kK ~ 1 , 

i.e. 

k(M x — rc T + gK~ 1 rc T K^ 1 ) + (M x — rc T + gK~ 1 rc T K~ 1 )k = 0, 
where we have made use of the relation kK = Kk. Using proposition 3.1 this yields 

M x = rc T — gK~ 1 rc T K~ 1 , (5.2) 

i.e. 

M x = k,M + M k . (5-3) 

if we use the Sylvester equation (3.1). In a similar way, for the time evolution of M, we have 

M t = 4(k 3 M + Mk 3 ). (5.4) 
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5.2 Evolution of 


With the evolution formulas (5.1)-(5.4), we can derive the evolution of We make use of 

the auxiliary vectors defined in (3.33), i.e. 


(I + M> (2i) = K l r, (5.5a) 

(I + iW>( 2i+1 ) = kK l r. (5.5b) 

By them are expressed as 

Sf 24 *') =c T K j u { - 2i \ (5.6a) 

s ( 2 i+i, 2 j) =c T K j u ( ~ 2i+1 \ (5.6b) 

jS r(2i,2j+i) = c r K j kuW, (5.6c) 

5 (2i+i,2j+i) =c T K j ku( 2i+l \ (5.6d) 


Taking x-derivative on (5.5a) we have 

M x u (2i) + (I + M)u {2i) = K l r x = K‘kr , (5.7) 

and further, by substitution of (5.2), we have 

(I + M)uf l) = K { kr - ( rc T - (5.8) 

which indicates the evolution of ■u ( ' 2l) in x-direction: 

U W = U (M+D - S^W 0 ) + 5 £( 2i ’- 2 V- 2 ). (5.9a) 

Similarly we can derive out the evolution of in x-direction 

t f+l) = U < 2i+2 > - ^( 2i +b0) u (0) + gS (2i + l,-2) u (-2) + 3eiW (2*) + ^( 2i -2) ; (5 . 9b) 




and the evolution of in t-direction: 

uf i] = 4 [W 2) (S (2 *’ 0) + 3 ei ^( 2i ’- 2 ) + 5 ^ 2i ’- 4 )) + 9 5( 2l ’- 2 V (0) + 3 ei u(" 2 ) + W" 4) ) 

- u^(S^ + 3 ei 5( 2i ’°) + 5 S( 2 *’- 2 )) + S^W 1 ) - 

+ 3eiu (2i+1) + gu^-V - 5 (2i ’ 0) (w (2) + 3ei-u (0) + gu^~ 2) )}, (5.10a) 

u? i+1) = %u(- 2 )(S( 2i+1 ’°) +3 ei ^ +1 -- 2 ) + gS^ i+1 ’~^) -g U (-Vs& +1 ’-V +U&+V 

- i i (°)(5 (2 *+ 1 ’ 2 ) + 3ei5 (2i+1 ’ 0) + gS {2i+1 ’~ 2) ) + (9e? + 2 5 )u (2i) + c/V 2i “ 4) 

+ 6e l5 u( 2i - 2 ) + 6 ei u( 2i+2 ) - S( 2i+1 ’°\uW + 3ei «(°) + W“ 2 )) 

+ 5 , ( 2 *+ 1 .i) tt (1 ) + gS^ 2i+l - 2 )( w (°) + 3eiu(" 2 ) + c/tt (—4) )]. (5.10b) 

These evolution of can be turned into the evolution of S^’^: 


g(2i,2j) = g(2i+l,2j) _|_ g(2i,2j+l) _ g(2i,0) g(0,2j) _|_ gg(2i,-2) g(-2,2j) ^ 

g(2i,2j+l) _g(2i,2j+2) _|_ g e ^g(2i,2j) _|_ gg(2i,2j-2) _ g(2i,0) g(0,2j+l) _|_ g(2i+l,2j+l) 

+ gS ^~2) s (-2,2 j+ l)^ 

s (2i+l,2j+l) = g(2i+l,2j+2) _|_ 3g 1( 5'( 2 *+l>2j) _|_ gg(2i+l,2j—2) _ g(2i+l,0) g(0,2j+l) + g(2i+2,2j+l) 

+ 5 5(2i+l,-2) 5 (-2,2j+l) + 3ei 5(2i,2j+l) + 5jS (2i-2,2j+l) ) (5.H C ) 


(5.11a) 

(5.11b) 
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and 


s (2i,2 j ) =4 [ gS {~2,2j)( S {2i,0) + 3eiS ( 2i-2) + gS^~^) + gS ^’" 2) {S {0 ’ 2j ^ + 3d S { ~ 2 ’ 2j) 

+ g s(~ 4 ’ 2j) ) - s( 0 ’ 2j \s {2i ' 2) + 3eiS (2i ’ 0) + gS^ 2i ’~ 2) ) + fiff 2 *- 1 ) 5 ( 1 . 25 ) 

_ ggf-b 2 ?)^ 2 *,-!) -f- g(2j+3,2j) _|_ 3 eij gf(2i+l,2j) _|_ g£i(2i—l,2j) _ g(2i,0) (g(2, 2j) 

+ 3ei5 (0 ’ 2i) + ff 5 ( - 2 ’ 2i) ) + S (2i ’ 2j+3) + 3eiS (2i ’ 2j+1) + s^*'" 1 )], (5.12a) 

sf ’ 2j+1) =A[gS^ +1 \S^ +3 ei 5( 2i >- 2 ) + 5 5( 2i -- 4 )) + ^ 2i ’- 2 ) (^ +1 ) 

+ 3 eiS ( " 2 ’ 2j+1) + fl 5 ( - 4 ’ 2j ' +1 >) - s(o, 2 j+i)(£( 2 i, 2 ) + 3 ei 5 ( 2 i,o) + g s^- 2 )) 

+ 5 ( 2 i> 1 ) S , ( 1 > 2 -?+ 1 ) - ^(-baj+i) S^ 2i “ 1} + S'( 2 *+ 3 . 2 J+ 1 ) + 3eiS' (2i+1 ’ 2j+1) 

+ 5<S , ( 2i_1,2j+1 ) — 5 , ( 2 *’°)( ( S'( 2 ’ 2 J+ 1 ) + 3eiS , ( 0,2: ' +1 - ) + ggl -2 ’ 2 -? 4-1 )) + 

+ 6ei5 (2i ’ 2j+2) + (9e? + 2 g)S {2i ’ 2j) + 6 ei gS^ 2i ’ 2j ~ 2) + g 2 S^ 2j ~% (5.12b) 

5f + l’2^ + l) =4[ & 5(- 2 >2j+l)(5(2 < +l-0) + 3 ei 5( 2 i +l.-2) + 5 S(2i+l.-4)) + 3 S (2i+lr2) (S (0,2j+l) 

+ 3ei5 ( " 2 ’ 2j+1) + 5 < S ( “ 4 ’ 2j+1) ) - S ( 0 ’ 2 j ' + 1 )(S (2i+1 ’ 2 ) + 3 ei 5 (2i+1 ’ 0) + g s ( 2i+1 ’~ 2) ) 

+ 5 (2»+l,l) s (l,2j+l) _ s (-l,2j+l) s {2i+l,-l) + 5 (2j+4,2j+l) + 6 eiS (2i+2 ’ 2j ' +1) 

+ Ge^ 2 ^ 1 ’ 2 ^ + (9e 2 + 2 g)S^ 2i+1 ^ + 6e l9 S^ +1 ’ 2 ^ + 5 2 s( 2 *+b 2 J- 4 ) 

— 5( 2 *+ 1 '°)(5( 2 > 2 i+ 1 ) _|_ 3eiS^ 0,2:,+1 ^ + gS^~ 2,2 ^ +l ^) + gf 2 ^ 1 - 2 ^ 4 ) 

+ (9e? + 2 5 )5 (2 * +1 ’ 2 ^ + 6 e l 5 S (2i+1 ’ 2j - 2) + g 2 S {2i+1 ’ 2j ~ 4) ]. (5.12c) 

It is not necessary to write out the S^ +1 ’ 2 ^ and gf 2 * 4-1 ’ 2 ?) due to the symmetry property 
= S^ ,l \ One can repeatedly use (5.11) and easily get higher-order x-derivatives of S'(O) 
by means of computer algebra, e.g. Mathematica. 

These derivatives of S^O) bring the following epKdV system [1] 


u t =u xxx + 6u 2 x - 6gs 2 x , (5.13a) 

$t —$xxx T 6'U j xSx QgSxhxi (5.13b) 

h t =h xxx + 6s 2 - 6 gh 2 x , (5.13c) 

Vt — v xxx T 6 v x u x T 6 gs x S x ’ \ (5.13d) 

w t =w xxx + 6s x S'). 0 ’ 1) - 6 gw x h x , (5.13e) 

where u,v,s,w,h are defined as in (4.5), and 

S { ~ 1 '- 2) = ^(h x +s 2 - gh 2 ), (5.14a) 

S(°’V = ^( Ux + u 2 - gs 2 ), (5.14b) 

which are from (5.11a). 


To finish the derivation of (5.13), one needs long and tedious verification in which the 
recurrence relations (3.30) are successively used. For example, for the first equation (5.13a), 


19 


after substitute the expressions of ut , u xxx ,u x ,s x one has 
g(^t ^xxx 6gs x ) 

= 5 2 5 (°,-2)[ 5 (0,- 4 ) __ ^-l,-2)5(0,-2) + 5 (-2,-2) ( _ 1 + £-(0,-1))] 
+ ff[ _ 5 (0,-2) ( _ 5 (0,0) (1 + £-(1,-2)) + £-(2,-2) + £,(0,-2)£.( 1 ,0) ) 
_ 5 (1.-2) + 5(0,-D(l + 5(h-2)) _ £.( 1 ,- 1 )£.(0,-2) ] 

_ 5(i,o) 2 + 5(0,0) 5(1,1) _ 5(2,1) + 5(3,0) _ 


It vanishes in light of the relations (3.30) with (i,j) = (0,-2) and (0,1). Equations (5.13b)- 
(5.13e) can rigorously be derived in a similar way but the procedures are long. Here we skip 
them. 


In [1] 

w 

A = —u + — (5.15) 

is introduced, by which the epKdV system (5.13) yielded the following coupled equation, (i.e. 

( 1 - 2 )) 


st =4 s xxx + 6 s x [i?(s 2 ) - A 2 — - —h£] ; (5.16a) 

s s 

6s 

At =4 A xxx — 6 A 2 A x + 6 A x R(s 2 ) - (R(s 2 )) x , (5.16b) 

s 


with the elliptic curve R(x) given in (1.3). This coupled system is integrable in the sense of 
admitting a continuous Lax pair [1] (also see Sec.5.3). 


5.3 Lax pair 

Let us consider (5.9) with i = 0. which reads 

u = u ^ — uu^ + gsu^~ 2 \ (5.17a) 

u ^ + gwu ( '~ 2 ' 1 + 3ei u^°\ (5.17b) 

After replacing u^ 2 \ u^ and S^LO) with (4.14), (4.16) and (5.14b), respectively, we have 

u^ = u 1 ' 1 ' 1 — uu^ + gsK~ 1 (wu — su^), (5.18a) 

= (K + gw 2 K~ l )u^ + (3ei — u x — u 2 + gs 2 )u ^ + {ul — gswK _1 )u^\ (5.18b) 


As in the discrete case, we can consider k to be its canonical form and then from the first rows 
of (5.18a) and (5.18b) we find the linear form 




-u+j^sw 

K + 3e + gs 2 — u 2 — u x + j^w 2 u — j^ws ) ’ 


(5.19a) 


where (j> is defined by (4.17). In a similar way from (5.10) we can find the time evolution of 0, 
which is formulated by 


4>t 


sl 0,1 ) u x \ , g ( (1 — vw)^ + v x w vs x — sv x 

Sx'^ 5i 0,1) K ( (1 - vw)^f - w(±=^) x -(1 -vw)^f-v x w 


(5.19b) 
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where 


Si 1 # = + uSi°’ 1] - gsw x , (5.20) 

and S^ 0,1 ) is given by (5.14b). 

(5.19) can be viewed as a Lax pair of the system (5.13), which can also be derived from the 
direct linearization approach [1], The compatibility gives equations (5.13a), (5.13b) and (5.13e). 


6 Straight continuum limits 

The skew continuum of the elpKdV system (4.7) was considered in [1] when they studied initial 
value problems of (4.7). Such a limit is performed by introducing the skew-change of variables 
(n, m) i—>■ (IN’ = n + m, m). 

Let us consider the straight continuum limit, where we first take 

777- 

m oo, 6 —oo, while — = r — To ~ 0 ( 1 ), ( 6 . 1 ) 

5 

with tq being a constant. Define 


^ — ^n,m —• Un(x'), $ — Sn,m —• S n (T), W — W n ,m —• ^n(x). 


( 6 . 2 ) 


Then, applying the Taylor expansions into (4.7) at r, the leading term (in terms of 1/6) of each 
equation yields the following semi-discrete equations 


d T (u + u) = 2 a(u — u) — (u — u) 2 + g(s — s) 2 , 

d T (ss) = (s — s) (as + as — w + w) + us 2 + us 2 — ss(u + u ), 

, W., _ W. 2 r,/ 

(a + u — — )(a — u -|- ) = a —Kiss), 

s s 

dr(u + -) + (u--) 2 = R(s 2 ), 
s s 

and we note both (4.7b) and (4.7c) yield (6.3b) in the continuum limit. 

For the full limit of (6.3), first we take 

n oo, a oo, while — = £ ~ 0(a 2 ) 

a 

and then introduce continuous variables x and t by 


x = t + £, 


t = 


s 

12 a 2 ’ 


(6.3a) 

(6.3b) 

(6.3c) 

(6.3d) 


(6.4) 


(6.5) 


with £ as an auxiliary variable. Then, under the coordinates ( x , t) both (6.3c) and (6.3d) yield 


+ +(u--^) = R(s 2 )-, ( 6 . 6 ) 

(6.3a) yields (5.13a); (6.3b) gives (5.13b) and for that we need to make use of a relation s xx = 

2 (gsh x — us x + w x ), which is obtained as a continuous limit of the summation of (4.6c) and 
(4.6d). 
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If we employ the transformation A = —u + ^ in the equations (6.3b) and (6.3c), it turns out 
to be 

d T (ss) = (s — s)(as + as — As + As) + us 2 + us 2 — sJ(u + u), (6.7a) 

(a + u — A — u)(a — u + A + u) = a 2 — R(slP). (6.7b) 

Their continuum limit gives the coupled system (5.16), where we need to use the relation u x = 
^( R(s 2 ) — A 2 — A x ), which is (6.6) in terms of A and u, derived from (6.3d). 

7 Conclusions 

In this paper a new class of solutions of the elliptic KdV systems (both the discrete (1.1) and 
the continuous (1.2)) has been uncovered. We made use of Sylvester-type equation with elliptic 

ingredient. Solutions can be classified by the canonical form of k, which are much richer than 

pure solitons. 

A Cauchy matrix dressed by dispersion relations usually satisfies a Sylvester equation. Start¬ 
ing from the Sylvester equation and dispersion relations, not only integrable equations can be 
derived but also their solutions and Lax pairs can be constructed. Such the Cauchy matrix 
approach is particularly powerful in the study of discrete integrable systems (see [6,8]), as well 
as continuous systems [7]. Dressed Cauchy matrices also play key roles in the so-called operator 
method [14-19], trace method [20], etc. 

Here we specially mention the connection with Inverse Scattering Transform (1ST), which 
is to express solutions of the Gelfand-Levitan-Marchenko equation via a triplet (A; b; c) where 
matrix A and vectors b and c satisfy some Sylvester equations [21-23]. Solutions are classified 
by the canonical form of A. The obtained breather-like solutions can describe potentials on 
half-line and with non-zero refection coefficients [21]. We believe such connections can be found 
for the fully discrete systems which were solved via 1ST recently [24,25]. Besides, since there is 
deep connection between Cauchy matrix approach and direct linearization scheme, it is possible 
to develop a more general direct linearization scheme. 

Adding elliptic information to integrable systems, either to equations or to solutions, is an 
interesting topic, which brings the study of the integrable systems into a larger area and more 
insight [26-33]. In this paper, we have shown that the Cauchy matrix approach works for the 
study of some elliptic integrable systems, i.e. some terms of these systems are formulated with an 
elliptic curve. The Sylvester equation (1.6) is our starting point. We derived the discrete as well 
as continuous elliptic KdV systems. Apart from solutions, we also obtained Lax representations. 
As for solutions, the discrete plane wave factor (B.2) and continuous one (B.6) are defined with 
the wave numbers ki which together with ATj obeys the elliptic curve (1.3). For the Lax pairs 
(4.18) and (5.19), (k, K ) plays the role of spectral parameters which also obeys the elliptic curve. 
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A Lower triangular Toeplitz matrices 


Here we collect some properties of Lower triangular Toeplitz (LTT) matrices. 
A JV-th order LTT matrix is a matrix of the following form 


r |wi ({a.}i ) = 


Let 


a\ 

0 

0 

0 

0 

\ 



02 

Oi 

0 

0 

0 




03 

02 

Ol 

0 

0 



(A.l) 

oat 

OAT-1 

OAT—2 ' ' ' 

02 

Ol 

/ 

NxN 


7 [Nl 

II 

iM?)h 





(A.2) 


then we have AB = BA, \/A,B £ T [JV] , i.e. T [JV| is a commutative set w.r.t. matrix multiplica¬ 
tion. Particularly, the subset 


Tf = {F £ T [wl | det(F) + 0} 


(A.3) 


is an Abelian group. 

Obviously, the Jordan block matrix 


r l f{a) = 


(a 0 0 
1 a 0 
0 1 a 


0 0 \ 
0 0 
0 0 

1 a ) 


(A.4) 


\ o o o • 

is a LTT matrix, and one can verify the following. 

Proposition A.l. If A £ Cj\r X N and T [ f\a)A = AT^^a), then there must be A £ 7 lN h 


If a,j £ C, then the LTT matrix (A.l) can be generated by certain function. Suppose that 
f{k) is an analytic function. Using Taylor coefficients 

Qj = tt 3 ^ dT 1 f{k)\ k =k 0 , 3 = 1,2, - - - ,iV (A.5) 

we can generate a LTT matrix. The matrix (A.l) with (A.5) is called a LTT matrix generated 
by f(k ) at k = ko, denoted by T lN] [f(ko)}, and f(k) is called generating function. In this sense, 
the Jordan block (A.4) is generated by f(k) = k at k = a, and the unit matrix I is generated 
by f{k) = 1. On the other hand, for any LTT matrix (A.l) with a,j £ C, it can be generated by 
the polynomial 

N 

a(k ) = ^2a>j(k — k^y~ l (A. 6 ) 

3 = 1 

with {oj} as coefficients. Next, by [/(&)][^ we denote a set of functions (equivalence class) in 
which all the functions have the same (N — l)-th order Taylor polynomial at k = ko as f(k) 
has. Say, f(k) ~ g(k) if they have same ( N — l)-th order Taylor polynomial at k = k$. Thus, 
the LTT matrix (A.l) can be generated by any f{k) £ [a(fe)]^. With such correspondence, we 
have the following. 
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Proposition A. 2. If A = T w [/(fco)] and B = T [N, \g(ko)\, then 

C = AB = T"[f(k 0 )g(k 0 )}. (A.7) 

i.e., AB is a LTT matrix generated by f{k)g(k) at k^. As a result, we have 

flT^[f j (ko)}=Tm[flf j (k 0 )}, 

j =i j =i 

and 

(T^lfiko )})- 1 = TM[l/f(ko)] 

iff(k o)/0. 

Proof. We only need to prove (A.7). Suppose 


A = (aij)NxN, B = (bij)NxN, C = (Cij)NxN- 
Then we have (with k = ho) 

A _ t(u\ ( _I_ fp-i. 


jpfj\ d k f(k), i>j , , bij = \ f^JY d k S(k), i>j 


0 , 


* < J , 


0 , 


i < J, 


and (with k = /cq) 


N 


C ij ~ 


( ai s b s j 


S=1 

N 


-77 d^gik) 


5(^™-( S -7)!' 

S J——T)u\ d k 3 l fW • 4s(*0> (* > i)> 


l=i-j 


1 ^ 

(-—yy I] J (f(k)g(k)), (i>j), 


l=i—j 


and Cjj = 0 when i < j. Thus, (A.7) is proved. 


□ 


In addition to the LTT matrices, we define the following skew triangular Toeplitz (STT) 
matrix: 

( b\ • • • &7V-2 fi/V—l bN \ 

b 2 • • • bN -1 0 

63 • • • 67 V 0 0 




\ t>N 


(A.8) 


0 0 0 / 


NxN 


The following property holds [7]. 
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Proposition A.3. Let 


j [N] = { H m({ bj }*)}. 


(A.9) 


Then we have 


(1) . H = H T , \/H G 7 !N] ; 

(2) . HA = ( HA) t = A t H. VA G 7 m , VH € T w . 


It can be extended to the following generic case. 


Proposition A.4. Let 



where 0 < Nj < N for j = 0, 1, • • • , N and Nj = N. Then we have 

(1) . AB = BA, VA, B € S w ; 

( 2 ) . h = h t , vh g gw- 

( 3 ) . HA = ( HA) t = A t H, VA € VH G S w . 


B Explicit forms of r and c 

Here we list out the explicit forms of r and c satisfying (4.1) and (5.1), respectively. 

B.l Solution to (4.1) 

(1). When r = r^({fcj}f) we have 


r = 


r ^({%}f) = (hd2, ■ ■ ■ An) T , with ri = pi, 


(B.l) 


where 



(B.2) 


and is a constant. 

(2). When T = r^fci), we have 



(B.3) 


where p\ is defined in (B.2). 
(3). When T = we have 


r l ? 2 ] (k N i+i) 
r l f 3 ] (k Nl+2 ) 


(B.4) 


V r l j Ns] (k Nl+{s _ 1} ) J 

where r\f l] ({fej}^ 1 ) and r [ j Ni] (kj) are defined as in (B.l) and (B.3), respectively. 
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B.2 Solution to (5.1) 

(1). When r = we have r = r [ ^({kj}^) is given in the form (B.l) and 

c = c l o({kj}i) = (ci,c 2) - • • ,c n ) t , with a = r h 


but here 

Pi = e^, & = kiX + 4 kft + with constant 

(2) . When T = r^Aq), we have r = r^(k\) is given in the form (B.l) and 

c = cjf^fci) = (ci,c 2 , • • • ,c n ) t , with a = r N -i, 
where p± is defined in (B.6). 

(3) . When T = r^ 1 we have r is given in the form of (B.4) and 

/ \ 

c l ? 2] 3 4 5 6 7 8 9 10 (k N i+i) 

c= c) JV3l (fc A r 1 + 2 ) 

\ c) Wsl (A: A r 1+(s _ 1) ) / 

where c [ „ Vll ({A.'j}) Yl ) and c [ f i] (kj) are dehned as in (B.5) and (B.7), respectively. 


(B.5) 


(B.6) 


(B.7) 


(B.8) 
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